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First-principles evaluation of Multi-valent cation insertion into Orthorhombic V 2 0 5 
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A systematic first principles evaluation of the insertion behavior of multi-valent 
cations in orthorhombic V2O5 is performed. Layer spacing, voltage, phase stability, 
and ion mobility are computed for Li + , Mg 2+ , Zn 2+ , Ca 2+ , and Al 3+ intercalation in 
the a and 5 polymorphs. 


A promising and realistic strategy to improve the en¬ 
ergy density beyond the capability of current Li-ion bat¬ 
tery technology is to transition to a battery architecture 
based on shuttling multi-valent (MV) ions (e.g. Mg 2+ or 
Ca 2+ ) between an intercalation cathode host and MV 
metal anodePf 2 Specifically, improvement in the volu¬ 
metric energy density arises from the combination of us¬ 
ing a multi-valent metal as the anode as opposed to an 
insertion structure (e.g. ^ 3833 mAh/cm 3 volumetric 
capacity for Mg metal compared to ^ 800 mAh/cm 3 
for graphite), and storing more charge per ion in the 
cathode P® 

One of the major bottlenecks preventing the develop¬ 
ment of MV battery technology, however, is the poor elec¬ 
trochemical performance of potential MV cathode mate¬ 
rials, thought to originate predominantly from poor MV 
ion mobility in the intercalation host structureP® More¬ 
over, the simultaneous challenge of developing function¬ 
ing MV anodes and electrolytes compatible with candi¬ 
date cathode materials has limited the ability to exper¬ 
imentally isolate and evaluate cathode electrochemical 
performance,^ and as such there is a general dearth of 
reliable data on MV ion intercalation in the literature to 
date to guide the ongoing search for new MV cathode 
materials with improved performance. 

Nevertheless, reversible electrochemical Mg 2+ interca¬ 
lation has been successfully demonstrated in a handful 
of cathode hosts, namely Chevrel MoeSg (~ 135 mAh/g 
capacity at ~ 1.0 — 1.3 V vs. Mg metal) P as well as 
layered V 2 O 5 (^ 150 mAh/g at ~ 2.3 — 2.6 V)P^ and 
M 0 O 3 (~ 220 mAh/g at ~ 1.7 — 2.8 V)PThe orthorhom¬ 
bic V 2 O 5 structure is especially interesting because it 
has also demonstrated the ability to reversibly intercalate 
Ca 2+ and Y 3+ in addition to Mg 2+ ions! 6 ^ First princi¬ 
ples calculations (described in more detail in the supple¬ 
mentary informatiorl 29 ^) have proven to be an accurate 
and effective method to systematically assess the elec¬ 
trochemical properties of Li-ion batteries ,®^ 0 and have 
also been used to study the process of ion intercalation 
in layered materials, such as graphite 11 and V 2 O 5! 12 ^ 



FIG. 1 : a) The V2O5 structure of both the a and 5 
polymorphs on the b — c plane with the yellow spheres 
indicating the intercalant sites while b) shows the a 
and 6 polymorphs on the a—b plane. As indicated by 
the dashed blue regions, both the polymorphs differ 
by a change in the stacking of the V2O5 layers. 


In this work, we have performed a systematic first prin¬ 
ciples study of MV ion intercalation in the orthorhombic 
a— and S— V 2 O 5 polymorphs by evaluating the struc¬ 
tural change, voltage, thermodynamic stability, and in¬ 
tercalant mobility for Li + , Mg 2+ , Zn 2+ , Ca 2+ , and Al 3+ 
insertion and comparing to data in the literature when 
available. 

The crystal structure and intercalation sites of the a— 
and 5— V 2 O 5 polymorphs 16 19 ^ are shown in Fig. [I] Per¬ 
pendicular to the 6 -axis (i.e. in the a—c plane), the or¬ 
thorhombic V 2 O 5 structure consists of layers of alternat¬ 
ing corner— and edge—sharing VO 5 pyramids (shown in 
red), each consisting of 4 V—O bonds that form the base 
and one short V=0 bond that forms the apex. The in¬ 
tercalation sites (yellow spheres) are situated in between 
the layers, and assuming no limitation in the number of 
redox centers, the theoretical gravimetric capacities for 
AV 2 O 5 where A = Li, Mg, Zn, Ca and A1 are 142, 260, 
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217, 242 and 385 mAh/g, respectively. Structurally, the 
main difference between the a and S polymorphs is a 
shift in the layer stacking, indicated by the dashed blue 
lines in Fig. [TJd, with alternate V 2 O 5 layers displaced in 
the a-direction by half a lattice spacing, accompanied by 
a change in the interlayer distance and the anion coor¬ 
dination environment of the intercalation sitesP^ While 
8 oxygen atoms coordinate the intercalant ion in a (for 
Mg, there are two Mg—O bonds with length 2.11 A, 
two with ^ 2.39 A, and four with ~ 2.46 A, respec¬ 
tively), “4+2” oxygen atoms coordinate the intercalant 
in S (for Mg, there are four Mg-0 bonds with length 
~ 2.05 — 2.07 A, and two with ~ 2.33 A). 

In Fig. UK the interlayer spacings in the a and S 
polymorphs (filled and hollow bars, respectively) are 
shown for empty V 2 O 5 and intercalated AV 2 O 5 , where 
A = Li, Mg, Zn, Ca, and Al. To better capture the 
increased effect of van der Waals effects in the deinter- 
calated limit, the interlayer spacings for empty V 2 O 5 
(4.46 A for 5.03 A fo r 5) are calculated using the 
vdW-DF2 functiona l 20 * 21 * rather than standard DFT as 
the latter significantly overestimates this spacing (4.75 A 
for a ; 5.2 7 A for (5) compared to experiment (4.37 A for 
a). I 12 * 15 * 1 - As detailed in the supplementary information, 
Al 3+ intercalation in the a— V 2 O 5 structure is found to 
be mechanically unstable and relaxes to the 5 polymorph 
in our calculations, and we therefore remove it from fur¬ 
ther consideration in this study. 

At the same intercalant composition, the S structures 
consistently have larger layer spacings than <a, ~ 3 — 5 % 
larger for Li, Mg, and Zn and ^ 10 — 12 % for Ca and 
empty V 2 O 5 . With the exception of Ca intercalation, 
which increases the layer spacing by more than 10 % 
in both polymorphs, the change in the layer spacing is 
much smaller in 5 than <a, less than 2 % for Li + , Mg 2+ , 
Zn 2+ , and Al 3+ intercalation in 5— V 2 O 5 compared to 
- 9 - 14 % for Li+, Mg 2 +, Zn 2 +, and Ca 2+ in a-V 2 0 5 . 
The behavior for Ca 2+ is consistent with intercalation 
in the spinel system, 2 ^ where the volume change is also 
much larger than for Li + , Mg 2+ , Zn 2+ , and Al 3+ inter¬ 
calation, and in general may be attributed to the larger 
ionic radius of Ca 2+ in comparison to the other ions . 23 
Al 3+ intercalation in 5— V 2 O 5 , in contrast to the other 
ions considered, is accompanied by a contraction of the 
layers, which is consistent with its small ionic radius and 
higher positive charge density that strengthens the at¬ 
traction with nearby oxygen ions. 

The average voltages of the compounds computed us¬ 
ing the method of Aydinol et a/P^are plotted in Fig. If 
and are referenced to the potential of the bulk metal of 
the corresponding intercalating ion (i.e., Li metal for Li + 
intercalation, etc.). The average voltages computed for 
Li, Mg, and Ca intercalation compare very well to avail¬ 
able experimental data: ^ 3.2 — 3.4 V for Li measured 
by Delmas et al.^l ~ 2.2 — 2.4 V for Mg measured by 
Gershinsky et aZ.P and ~ 2.4 — 3.1 V for Ca measured 
by Amatucci et a/P In general, the Li polymorphs have 
the highest voltage, followed by Ca, Mg, Al, and Zn, 



FIG. 2 : a) Plots the layer spacing values for the empty 
and intercalated versions of AV2O5 (A = Li, Mg, Ca, 
Zn and Al) for both the a and S polymorphs, b) 
Displays the calculated average voltage values for the 
intercalation of the different ions and c) shows the 
energy above hull, which quantifies the stability of a 
structure, for the empty and intercalated versions of a 
and S. The filled regions in all the graphs correspond 
to the a structure while the hollow regions correspond 
to the 6 structure. Note that the energy above hull 
for o-CaV 2 0 5 is 0 meV/atom, implying that it is a 
ground state configuration in the Ca-V-O system. 
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which reflects both the same order and approximately 
the same potential difference indicated by the electro¬ 
chemical series (—3.04 V vs. SHE for Li, —2.86 V for 
Ca, —2.37 V for Mg, —1.66 V for Al, and —0.76 V for 
Zn). In comparison, the voltage difference between the 
V 2 O 5 polymorphs is much smaller for a given intercala¬ 
tion chemistry. For Li, Mg, and Zn the insertion voltage 
is higher in S (3.36 V, 2.56 V, and 1.09 V, respectively) 
than in a (3.18 V, 2.21 V, and 0.68 V), unlike for Ca 
where a is higher (3.13 V for a; 3.02 V for (5). 

Fig. [ 2 ): displays the energy above the convex ground 
state energy hull (E A hull) of the deintercalated and inter¬ 
calated V 2 O 5 polymorphs with respect to the intercalant- 
V-0 ternary phase diagram. The ternary ground state 
hulls were determined from the available calculated com¬ 
pounds in the Materials Project database.®! A predicted 
thermodynamically stable structure will have a E A hull 
value of 0 meV/atom while higher (more positive) E A hull 
values indicate greater instability, which may be reflected 
in experimental difficulties in synthesis or decomposition 
during battery operation. Note that the E A hull values 
calculated here reflect the ground state (i.e. 0 K), and 
entropy contributions, which scale with k^T, can stabi¬ 
lize certain structures at higher temperatures. 

In the deintercalated limit, V 2 O 5 is thermodynami¬ 
cally stable in the a phase, but S is only ^ 13 meV/atom 
higher in energy, indicating the possibility of metasta¬ 
bility at room temperature. For Li intercalation, the a 
and S structures are 82 meV/atom and 57 meV/atom 
more unstable than the ground state orthorhombic 
7 —LiV 205 structure, which has a different orienta¬ 
tion of the VO 5 pyramids®! along the c-direction 
shown in Fig [jji, but the S structure can remain 
metastable and has shown to be reversibly cycled 
electrochemically.®! S— MgA^Os, which has been synthe¬ 
sized experimentally,® is only ^27 meV/atom more un¬ 
stable (compared to ~ 102 meV/atom for a) than the 
thermodynamic ground state, a two-phase equilibrium 
consisting of MgVOs and VO 2 . Similarly ci-ZnV^CF is 
only ^31 meV/atom more unstable than the ground 
state (ZnO and VO 2 ), indicating that a metastable syn¬ 
thesis comparable to the Mg system may be possible. As 
Al intercalated a+V^Os displays mechanical instability 
in our calculations, when relaxed its energy is not de¬ 
fined, but the Al intercalated ci-phase is ~ 158 meV/atom 
unstable compared to the ground state ternary equilib¬ 
rium of AI 2 O 3 , VO 2 and V 3 O 5 . With the exception of 
a— CaV 205 , which is the ground state in the intercalated 
Ca-V 205 system, the S structures tend to be more stable 
than a in the discharged state (by 25 meV/atom for Li; 
75 meV/ atom for Mg; and 91 meV/atom for Zn), and 
accordingly the insertion voltages for 5 are higher than 
a for Li, Mg, and Zn insertion but lower for Ca inser¬ 
tion, as observed in Fig [2}d. Given that the intercalant 
sites in a and S are coordinated by 8 and “4+2” oxygen 
atoms respectively, the stability of the discharged (LV 2 O 5 
structures for Li, Mg and Zn, and ce-A^Os for Ca align 
well with the preferred coordination environment of the 


Path Distance (%) 



FIG. 3 : The activation barriers for the diffusion of the 
different intercalating ions in the a and 5 polymorphs 
are plotted in a) and b) respectively. The solid lines 
correspond to the empty lattice limit (charged state) 
while the hollow lines correspond to the full lattice 
limit (discharged state). 


respective ions, as tabulated by Brown. 26 !Hence for inter- 
calant ions that prefer a lower coordination number (i.e., 
coordinated by a maximum of 6 neighboring atoms), an 
a -+ S transition upon insertion in V 2 O 5 is likely. 

Fig. I displays the migration energies for intercalant 
diffusion along the a-direction in the a (Fig. and 
5 (Fig. [ 3 ) 3 ) polymorphs plotted against the normalized 
path distance calculated with the Nudged Elastic Band 
method. 27 The solid lines correspond to migration ener¬ 
gies obtained in the empty lattice limit (charged state), 
and the dashed lines correspond to the fully intercalated 
limit (discharged state). As elaborated upon in the sup¬ 
plementary information, converging the migration ener¬ 
gies in structures that exhibit a high degree of thermo¬ 
dynamic instability may not be possible, as was the case 
for Li, Mg, and Zn in the intercalated a— V 2 O 5 struc¬ 
ture, and for Ca in the intercalated 5— V 2 O 5 structure. 
In lieu of determining the Mg migration barrier in the 
fully discharged a— V 2 O 5 structure, we have computed 
the energy for Mg migration in a half intercalated struc¬ 
ture with a specific ordering of Mg ions, referred to as the 
“e” phase, which has also been observed in the Li-A^Os 
system . 16 

In Fig [3j the maximum energy difference encoun¬ 
tered along the diffusion path defines the migration bar¬ 
rier (E m ), which provides an approximate estimate of 
the ionic diffusivity. As a guide, at room tempera¬ 
ture, E m ^ 525 meV corresponds to a diffusivity of 
^ 10 -12 cm 2 s -1 , and a 60 meV increase (decrease) in 
the migration energy corresponds to an order of magni¬ 
tude decrease (increase) in diffusivity. Due to stronger 
interactions between a multivalent intercalant and the 
surrounding anion environment, the migration barriers 
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within the same host structure, for example Al 3+ , are 
generally higher than the divalent ion barriers (Mg 2+ , 
Zn 2+ , Ca 2+ ), which are generally higher than the barri¬ 
ers for Li + . For the divalent intercalants, the trend in the 
migration barriers is Ca 2+ 1700—1900 meV) > Mg 2+ 
975—1100 meV) > Zn 2+ 305 meV) in the 
aephase, but Mg 2+ 600—800 meV) > Zn 2+ 
375—425 meV) > Ca 2+ 200 meV) in the S 
phase. The energy above the hull (Fig [2]:) ranked 
from the lowest to highest reflects this same trend, with 
Ca 2+ > Mg 2+ > Zn 2+ in a and Mg 2+ > Zn 2+ > Ca 2+ for 
6, and highlights the positive correlation between high in- 
tercalant mobility and low thermodynamic stability. For 
both V 2 O 5 polymorphs considered, the change in the mi¬ 
gration barrier from the deintercalated to intercalated 
limit for the same diffusing species is much smaller than 
the variation across intercalating ions. 

Although the a and S polymorphs of V 2 O 5 are struc¬ 
turally very similar as earlier discussed, the anion co¬ 
ordination environment and therefore diffusion topology 
of the migrating intercalant vary significantly, which ac¬ 
counts for the different shape of the migration energies 
seen in Fig and Fig [3 ]d. In the a phase, the stable 
insertion site is coordinated by 8 oxygen anions which is 
connected to the adjacent insertion site along the a-axis 
by a 3-coordinated shared face. The shape of the mi¬ 
gration energies shown in Fig. [3^i, therefore, reflect the 
change in coordination of 8 —>3 —)>8 encountered by the 
diffusing species with the migration barrier correspond¬ 
ing to passing through the shared face. For the S phase, 
the stable insertion site adopts a “4+2” coordination and 
shares a corner with the adjacent insertion site along the 
a-axis. To migrate to this site, the intercalant passes 
through a 3-coordinated face shared with an intermedi¬ 
ate 5-coordinated (pyramidal) site, and finally performs 
a symmetric hop to the next insertion site. The change 
in the anion coordination along the diffusion path is then 
“4+2” -^3-^5—>3-^ “4+2”, where occupation of the inter¬ 
mediate pyramidal site corresponds to a local minimum 
in the migration energy, as is reflected in Fig[ 3 ] 3 . Overall, 
the migration barriers are also lower in the S phase com¬ 
pared to a (significantly lower for some cases), which 
we attribute in large part to the smaller coordination 
change during the migration process encountered in S. 
Also, the change in the relative order of the migration 
barriers of divalent ions between a (Ca > Mg, Zn) and S 
(Mg, Zn > Ca) can be explained by the correlation be¬ 
tween the “preferred” coordination environments of the 
respective ions and the available anion coordination en¬ 


vironments around the intercalation sites.® In a given 
structure, migration barriers are higher for an ion whose 
preferred coordination aligns with that of the coordina¬ 
tion environment available for the intercalant site com¬ 
pared to an ion whose preferred coordination is different 
from that present in the structure. For example, Ca is in 
its preferred 8 -coordinated site in a and hence has higher 
barriers than Mg and Zn, which are not in their respec¬ 
tively preferred 6 and 4 coordinated sites. Whereas in S, 
Ca is present in an unfavored “4+2” coordinated site and 
hence has lower barriers than either of Mg or Zn, which 
are closer to their preferred coordination environments. 
Our results thus lend support to the hypothesis that co¬ 
ordination of the intercalation site is a good screening 
criterion for identifying fast multi-valent cation diffusers. 

An ideal MV cathode intercalation host must pos¬ 
sess several properties—high capacity, high insertion volt¬ 
age, and MV ion mobility, while simultaneously minimal 
structural change and thermodynamic instability. From 
the systematic first-principles study performed in this 
work, we are able to evaluate all of the candidate ma¬ 
terials across each of these criteria. On the basis of ion 
mobility, Al 3+ intercalation appears unfeasible at room 
temperature in V 2 O 5 due to its prohibitively high migra¬ 
tion barriers, and although Zn 2+ intercalation is deter¬ 
mined to be facile in both polymorphs and relatively sta¬ 
ble in the S phase, the insertion voltage is low. Mobility of 
Mg 2+ and Ca 2+ is determined to be poor in the a phase, 
but intercalation of these ions in the S phase appear most 
promising, with sufficiently high voltage (3.02 V for Ca, 
and 2.56 V for Mg) and mobility (E m ~ 200 meV for 
Ca and ~ 600—800 meV for Mg) albeit with moderate 
thermodynamic instability (27 meV/atom for Mg and 
40 meV/atom for Ca above the ground state hull in the 
discharged state). 
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